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Abstract We study a model of particle acceleration coupled with an MHD 
model of magnetic reconnection in unstable twisted coronal loops. The kink 
instability leads to the formation of helical currents with strong parallel electric 
helds resulting in electron acceleration. The motion of electrons in the electric 
and magnetic fields of the reconnecting loop is investigated using a test-particle 
approach taking into account collisional scattering. We discuss the effects of 
Coulomb collisions and magnetic convergence near loop footpoints on the spatial 
distribution and energy spectra of high-energy electron populations and possible 
implications on the hard X-ray emission in solar flares. 
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1. Introduction 


Most models of particle acceleration assume that particles gain energy in the 
corona and then are transported into the chromosphere where they produce non- 
thermal radiation. This assumption follows from the standard solar flare model 
(e.g. LShibatal 119961 [Yokoyama and Shibata| 1199811 . with the primary energy re¬ 
lease site located above the top of a flaring loop or loop arcade. However, the 
standard scenario faces a number of difficulties ( MacKinnon and Brown, 1989| 
|Benka and Holman, 1994| [Brown et al., 200^ . The main problem is the large 
flux of electrons required to produce the observed hard X-ray (HXR) intensity. 
Also, unless precipitating electrons are accompanied by ions, they may result in 
a strong return current, creating additional electrodynamic challenges, such as 
preventing the electron beam from reaching the chromosphere. 
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It has been suggested recently, that re-acceleration of electrons in the chro¬ 
mosphere may help to overcome the “number problem” dBrown et al, 2009| ). 

In the proposed scenario electrons are initially accelerated in the corona and 
transported to the chromosphere. However, unlike in the standard thick target 
model (see IBrownl II97I|1 . electrons lose their energy much slower due to re¬ 
energization by the electric field in the chromosphere. Indeed, this is possible 
even with comparably small field strength. In the presence of Coulomb colli¬ 
sions, the value of electric field E required to accelerate electrons depends on 
their energy Sg as E = qnl(eSe), where q is an almost constant parameter 
{q ~ O.OITcq (see e.g. |Spitzer[ [I9h2|l . Thus, for the typical chromo¬ 

spheric density of 10^^ cm“^, one needs a field of at least 30 V m“^ to accelerate 
electrons from 10 eV (approximately thermal energies), but only 10“^ V m“^ to 
accelerate electrons from 40 keV. This electric field can be achieved even with 
the classical resistivity provided the current density is about I A 

Particle acceleration in fragmented electric fields distributed inside flaring 
loops can provide another opportunity to solve the “number problem”, yield¬ 
ing acceleration efficiencies broadly consistent with those deduced from HXR 
observations. Such an acceleration model would help to reduce energy lost dur¬ 
ing particle transport and to avoid some fundamental electrodynamic issues 
characteristic to acceleration models involving a single localized current layer, in¬ 
cluding the return current effect (e.ff. [Knight and Sturrock|ll977|l or ion-electron 
separation (see e.g. [Zharkova and Gordovskyy 120041) . 

There are a few previous studies of distributed particle acceleration. [Vlahos, Isliker, and Lepreti[ 
(I2004P considered acceleration of particles by a system of fragmented unstable 
current sheets derived using the cellular automata model, with a particle motion 
described by continuous-time random walks. Although the motion of individual 
particles is oversimplifield, this model gives a good description of stochastic 
acceleration of particle population by a complex, fractal network of acceler¬ 
ation centres. ITurkmani et al.\ (|2006l) investigated particle acceleration using 
a stationary model of the stressed coronal magnetic field. [Bian and Browning[ 

(|2008p used the kinetic approach to study stochastic acceleration by fluctuating 
fragmented electric fields in a turbulent reconnecting plasma and determined the 
relation between the distribution of current sheets and resulting particle energy 
spectra. More recently, Gordovskyy and Browning (2011, 2012) considered pro¬ 
ton and electron acceleration in twisted coronal loops using a time-dependent 
three-dimensional field configuration (see also [Browning et aL[l20111l . These con¬ 
figurations are attractive as they are expected to be ubiquitous in the corona, as 
loops can be twisted by photospheric motions (see e.o. lBrown et a/. 1 120031) or can 
emerge already twisted ( [Hood, Archontis, and MacTaggart, 2012^ . In any case, 
they contain excess magnetic energy which can be released if a loop becomes 
unstable. 

Magnetic field convergence is expected near chromospheric footpoints of coro¬ 
nal loops and should have significant effect on particle motion. In the present 
paper we examine the effects of both magnetic field convergence and Coulomb 
collisions on particle acceleration and re-acceleration in twisted loops, effects 
which have not been considered in our previous works. 
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Figure 1. Initial magnetic field for the model with convergence near footpoint boundaries 
(panel a) and corresponding magnetic induction on the r = 0 axis (panel b). Panel (c) shows 
the final twist profile in the flux tube. 


2. Evolution of twisted magnetic loop and particle motion 


We consider twisted magnetic flux tube models with and without field conver¬ 
gence near footpoint boundaries. 

Magnetic reconnection in initially cylindrical flux tubes has been extensively 
studied using analytical initial equilibria ( [Browning et al., 2008[ Hood, Browning, and van der Linden, 2009 ) 
These models have been already employed to study particle acceleration in 
[Gordovskyy and Browning (12011112()12|l , so that we do not discuss them further 
here. 

A configuration with twisted magnetic field converging near footpoints can 
not be set up analytically. Therefore, it is derived numerically, using ideal MHD, 
by twisting the footpoints of an initially current-free magnetic flux tube. We use 
the simulation domain x = [—SLoiSLo]; U = [—SLoiSTo]; z = [—IOLq; lOTo]- 
In this model scaling parameters are the same as in [Gordovskyy and Browningj 


(|2012p : the scale length is Lq = 10® m, the characteristic magnetic induction is 


10 - 


,-3. 


kg m 
and the 


i3o=4xlO“®T and the characteristic density is po = 3.3 x 
which yields the characteristic Alfven velocity uq = 1.95 x 10® m s 
characteristic timescale to = 0-53 s. 

Initially, the field is potential, formed by two point ’’sources” located at ri = 
{0; 0; —10.5Lo} and r 2 = {0;0;10.5Lo} outside the simulation domain (Figure 
la): 


B(r) = Bs 


r - ri 

|r - ri|3 


r-r2 

|r - r 2 |® 


( 1 ) 


This yields a flux tube with the magnetic field at footpoints nearly 10 times 
higher than the field at the centre (Figure lb). The initial density and pressure in 
the MHD simulations are uniform: p{t = 0) = po and p{t = 0) = 5 x 10~^Pq^Bq, 
i.e. it is magnetically dominated plasma. 

The twist is introduced using the azimuthal velocity at the 2 = ±10To 
boundaries (corresponding to the photospheric footpoints). The angular velocity 
depends on the radius as a; ^ r^/rg exp(—r/rg), which is generally consistent 
with observed sunspot rotation (e.ff. [Gopasyuk[ [T977I1 . where rg = 0.2Lg. The 
corresponding twist (defined as t) = uj(r, t)dt) is shown in Figure 2c. The 
resulting configuration (at t = lOOto) has a maximum twist of approximately Ttt 
and is nearly force-free but is, apparently, unstable: the kink instability develops 
at t Ri 120to- 
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Figure 2. Evolution of magnetic field (green lines) and current density in the twisted flux 
tube with convergence near footpoints. Yellow and pink colours show \j\ = 4jo and \j\ = Sjo 
iso-surfaces, respectively. Times shown are in the units of £o. 




Figure 3. Variation of magnetic, kinetic and internal energies in twisted magnetic flux tube 
with the convergence near footpoint shown in Figure 2. Verical dash at t ^ lOOfo denotes the 
time when footpoint rotation stops. Energy is measured in units of BqL^. 
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Electron acceleration in twisted coronal loops 


It should be noted, that this time is much shorter than realistic timescale 
of development of kink instability. Indeed, twisted coronal loops can remain in 
equilibrium for up to 10^ — 10^ Alfven times. However, since we are interested in 
the magnetic reconnection occuring after the loop kinks, the development of the 
instability is speeded up by the velocity noise (see [Gordo vskyy and Browning! 

[MU)- 


Further evolution of the twisted flux tube once it becomes unstable is studied 
using resistive MHD simulations with the LareSD code ( |Arber et ai, 2001] ) (see 
also Gordovskyy and Browning pOlllI2012|l for more details). The resistivity in 
this part of the simulation is non-uniform and depends on local current density 
as 


v{j) 


0 j ^ jcr 
m j > jcr 


( 2 ) 


The critical current in the present simulations is jcr = and the 

resistivity value is rji = 10~^LQtQ^ (which is, essentialy, the inverse Lundquist 
number). 

Magnetic field and current density evolution for the model with magnetic 
convergence near footpoints is shown in Figure 2, while Figure 3 shows the 
variation of magnetic energy Em, kinetic energy Ex and internal energy Ej 
defined by following integrals over the domain volume: 


EmW = f B^dn 
2/ro J 

volume 

Exit) = \ J pv'^dn 

volume 

Blit) = J pdn, 

volume 

where dO is the elementary volume. The field and current for initially cylin¬ 
drical flux tube can be seen in Gordovskyy and Browning (2011, Figures 3 
and 4). It can be seen that both configurations (with and without magnetic 
convergence) demonstrate similar behaviour. In both cases, when the twisted 
flux tube kinks, the current density rapidly increases forming a helically-shaped 
structure. The decrease in magnetic energy at this stage is due to the dissipation 
of the azymuthal component of magnetic field and reconnection between twisted 
loop’s field lines and ambient field lines. Gradually, the current distribution 
becomes more filamentary and the current density decreases. Reconnection stops 
after about lOOto after kink instability occurs. During the rapid reconnection 
{t = 130 — 180 to, see Figure 3), the electric field in the system is order of 
10^ V m“^, although it occupies only a tiny fraction of the model volume (see 
Gordovskyy and Browning, 2012). The main difference between the initially 
cylindrical configuration and one with converging field is that, in the latter case, 
stronger magnetic field results in substantially enhanced current density near 
the footpoints. 
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Figure 4. Density profile within the flux tube used to account for the effect of Coloumb 
collisions. The “coronal” part (| 2 :| < 7.5 Mm) of the flux tube has a constant density of 
10® cm~® while the “chromospheric” part (|z| = 8 — 10 Mm) has an exponential denisty 
distribution with the scale height of about 170 km. 
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The obtained time-dependent electric and magnetic field configurations are 
used to calculate test-particle trajectories. Initially, protons and electrons in 
the model are uniformly distributed in the domain volume, with Maxwellian 
velocity distribution corresponding to T = 0.8 MK and an isotropic pitch- 
angle distribution. Since the scale-length of our MHD model is much larger 
than particle gyroradii, we use the relativistic guiding centre approximation to 
calculate particle trajectories (see [Northrop I1963|l : 


dr 

dt 

u 

d{iv\\) 

dt 

dfi 

dt 

7 


u -I- 




-E-b--(b-V(KB))-t- 
m 7 


V 


q 

6a 
St 




2/1 


pa 

2/1 an f Sa\ 


( 7 n||) an 
V V2 


- coll 


B(1 - (7n||)VV2) vv^ystj 


J coll 


I + (7f'||)^ + 2/iS 


(3) 

(4) 

(5) 

( 6 ) 
(7) 


where U|| is the particle velocity along the magnetic field (yfn is treated as a 
single variable for the sake of convenience), and u is the perpendicular velocity. 
The drift velocity u^; is defined as u^; = Also, b is the magnetic field 

direction vector b = B/S, /i is the magnetic moment defined as /i = 

7 is the Lorenz factor (7 = (1 — j k = \/l — u%lc^. 

The terms in [ ]coU brackets are to account for collisional energy losses and 
pitch-angle scattering, which are not usually incorporated in test-particle models. 
For typical coronal and chromospheric densities collisions are predominantly im¬ 
portant for particles with enegies in 1 — 100 keV range and, therefore, relativistic 
effects in the collisional terms can be ignored. 
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The experiments including Coulomb collisions are performed with the den¬ 
sity profile shown in Figure 4, where the region with increasing density (|z| = 
8 ... 10 Mm) represents the chromosphere. Collisional energy losses are taken into 
account through the decrease of the particle velocity V as follows: 


dV n 


( 8 ) 


where the particle velocity for the purpose of Coulomb collisions is defined as 
V = -I- Mg. (The drift velocity u is defined solely by local field parameters 

and represents predominantly bulk plasma flow, and, therefore, is excluded from 
the collisional terms.) Here n is the ambient plasma density and a is the constant 
a = 2q/rril; for typical chromospheric parameters q = 27re^ In A is approxi¬ 
mately 2.87 X 10“^^ eV^ cm^ ( Syrovatskii and Shmeleva, 1972[ ). Essentially, this 
parameter determines how deep an accelerated particle can penetrate into dense 
plasma; in case of an uniform density the stopping depth can be evaluated as 
follows: 

2 


f ^max \ 


1.75 


/ dd^nitial \ 


( ” V 

V 10®CTO“^ / 


\lMmJ \ IkeV J V10®cto“^/ 

Pitch-angle scattering is implemented through random changes of test-particle 
pitch-angle ratio a = cos6 = v\\/V as 


da 

dt 


= Ii{ dt 


1 — (a -I- Aa/2)' 


'V3 




Aa — liidt 


1 — (a — /S.a/2y 




Aq;2 


Aa, 


(9) 

where n(a;) is the operator yielding 1 or 0 with the probability x ov 1 — x, re¬ 
spectively. Provided Aa is sufficiently small, the pitch-angle distribution evolves 
as 


dN{E,a) 




^dN[E,a) 
da 


as expected. 

In all the considered test-particle simulations we use open boundaries (be. 
particles are free to leave the domain; their trajectories are calculated until they 
cross a boundary). 


3. Non-thermal electrons in a twisted coronal loop 

Here we mostly focus on the spatial distribution of non-thermal electrons, which 
have direct implications for observations. The energy spectra of particles accel¬ 
erated in initially cylindrical twisted magnetic flux tubes have been discussed 
by [Gordovskyy and Browning (I2()ll|l and [Browning et al.\ (I2()ll[) . In the present 
paper we compare those with the energy spectra of electrons accelerated in the 
presence of magnetic convergence and Coulomb collisions. 

The electron energy distributions (Figure 5) at the early stages of recon¬ 
nection are rather similar to those obtained in simulations with no magnetic 
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Figure 5. Energy spectra of electrons in the model with magnetic convergence and Coulomb 
collisions. Times (normalized by £o) shown in the panels. The shape of the distribution 
below 1 keV is not realistic on panels (b-f) (see text for details). 


convergence and no collisions (Gordovskyy and Browning, 2011). Thus, during 
the fast energy release stage {t = 130 — 150to) the spectra are combinations 
of a Maxwellian thermal distribution and nearly power-law high-energy tail. 
However, at the later stages {t = 150 — 250to), when electric fields are gradually 
decaying and the collisions become dominant, the part of the spectra around 
few keV becomes harder and, at some point, a gap appears between the thermal 
part and high-energy part of the spectra. This is similar to the energy spectra 
of electrons precipitating in the thick-target model. 

It should be noted, that the Equation (8) is valid only for non-thermal elec¬ 
trons {i.e. when £e ^ kT and, therefore, the shapes of spectra in Figure 5(b-f) 
below 1 keV are not realistic, although the total number of “thermal” electrons 
is correct. This should not affect the behaviour of particles at higher energies. 

Now let us consider spatial distribution of energetic electrons. The distribu¬ 
tions of electrons with energies > 5 keV along the flaring loop are shown in 
the Figure 6 for three different cases. They are calculated assuming the electron 
populations are symmetric in respect of the loop mid-plane z = 0. All three 
distributions correspond to the stage of fast energy release in the MHD model. 
It can be seen that in the case of initially cylindrical flux tube, high-energy 
particles are nearly uniformly distributed along the reconnecting loop, which 
can be interpreted as the result of rather uniform distribution of strong electric 
field “islands” along the loop’s length. The total number of particles in this case 
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Figure 6. Distribution of electrons with energies £ > 5 keV along the flux tube. Here z = 0 
corresponds to the centre of the flux tube and \z\ = 10 Lq corresponds to the footpoints. Panel 
(a) is for the initially cylindrical flux tube, panel (b) is for the case with convergence near 
footpoints, panel (c) is for the case with magnetic convergence and Coulomb collisions. All the 
distributions correspond to the stage of rapid energy release in the MHD simulations. 



Figure 7. Variation of cross-section of the volume occupied by high-energy electrons. Panel 
(a) is for the initially cylindrical flux tube, panel (b) is for the case with convergence near 
footpoints, panel (c) is for the case with magnetic convergence and Coulomb collisions. Solid 
lines are for the region \z\ = 8Lo"-10Loi dashed lines are for \z\ = 4Lo...6Loi dot-dashed lines 
are for central region | 2 | < 2Lo. Horizontal axis shows the time after onset of reconnection. (It 
should be noted that speed of reconnection in the models in the case (a) and in the case (b-c) 
can be slightly different and, therefore, times shown in panel (a) may correspond to slightly 
different phase compared to panels (b) and (c).) The shape of cross-section is determined 
visually (as it has quite sharp edges, see Gordovskyy and Browning, 2011), the size is defined 
as average diameter from two measurements, in x- and y- directions. 


is quite low because most electrons quickly get to one of the boundaries and 
leave the domain. 

In the case with converging magnetic field near footpoint boundaries high- 
energy electrons are rather uniformly distributed around the centre of the flux 
tube, but with some deficit close to the footpoints. This deficit is quite surprising, 
as one would expect a higher number of particles in these regions: particles are 
slowing down in stronger field and would spend more time near these boundaries. 
Apparently, the deficit is caused by more efficient acceleration due to stronger 
electric fields in these regions (see Section 2). 

In the case with magnetic convergence and Coulomb collisions the majority of 
the electrons are located near the footpoints, which indicates that test-particles 
spend most of their time moving through the dense plasma representing the 
chromosphere. This phenomenon might be explained as the result of simultane¬ 
ous action of three different effects: the converging magnetic field leads to high 
pitch-angles {i.e. Ug > un), which, along with collisional deceleration, prevents 
majority of electrons from escaping into central region of the flux tube or outside. 
At the same time, electrons do not thermalize completely as they are accelerated 
by direct electric field. 
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In our previous paper ( [Gordovskyy and Browning, 201^ we show that the 
cross-section of the volume occupied by accelerated particles increases with 
time. This happens due to the reconnection between magnetic field of twisted 
loop and the ambient magnetic field. Let us now examine how the cross-section 
sizes change with time at different distances from loop footpoints , in the three 
different models (Figure 7). 

It can be seen that in the case of an initially cylindrical flux tube the centre 
and the footpoints regions expand in the same way. Cross-section diameters 
increase from « ILq to ~ 2.5Lo within I20to- In the case with field convergence, 
the central region expands nearly at the same rate as the footpoint regions. Thus, 
around z = 0 the cross-section of the volume occupied by high-energy particles 
increases from w 1.2Lo to Ri 2.3Lo within about lOOto, while the region between 
l^l = 8...IOL0 expands from 0.4Lo to 0.7Lo- 

With Coulomb collisions taken into account, the footpoint sizes expand faster: 
the cross-section of the volume between |z| = 8...IOL0 increases from O.SLq to 
~ 1.2Lo compared with the increase from I.2Lo to 2.2Lq near the centre of the 
flux tube. One of the possible explanations for this effect is strong pitch-angle 
scattering leading to a wider spread of particles in radial direction. 

The radial expansion of the volume occupied by particles should result in the 
increase of the radial size of non-thermal hard X-ray sources observed in solar 
flares. Recent observations bv IKontar et Qi.l(l2010|l and |Kontar, Hannah, and Bian[ 
(|201ip show that this expansion may be happening in solar flares and the 
observed pace of expansion is comparable to what is seen in the numerical 
experiments. 


4. Summary 

In this paper we present preliminary results concerning electron acceleration 
in a kink unstable twisted coronal loop. For this purpose a model of twisted 
magnetic flux tube with the field convergence near its footpoints was developed 
using ideal MHD simulations. Further development of the kink instability re¬ 
sulting in magnetic reconnection and fast energy release was investigate using 
resistive MHD simulations with non-uniform resistivity. We formulated a set of 
relativistic gyro-kinetic equations taking into account energy losses and pitch- 
angle scattering due to Coulomb collisions and, based on the developed MHD 
model, studied electron acceleration during magnetic reconnection in twisted 
magnetic flux tube. 

The considered reconnection scenario appears to be an effective particle accel¬ 
erator. Taking into account scaling parameters (fo = 1/2 Tq ss 6 x 10^® J, 

to ~ 0.53s) once can estimate the total energy released during the reconnection. 
In the case of twisted magnetic flux tube with convergence near footpoints there 
is ~ 3 X 10^° J released during ss 50 s. The efficiency of particle acceleration, 
which strongly depends on the resistivity in these models, is quite low: approxi¬ 
mately 5—7% of the total energy released is carried by accelerated particles. This, 
along with the fact that the magnetic reconnection occurs in rather simple single¬ 
loop configuration, makes this model comparable to small self-contained flares 
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dAschwanden et ai, 2009| . We emphasize, however, that all these estimations 
are subject to the scaling parameters. 

The considered electron acceleration model is quite attractive because the 
particle acceleration is not localized at the top of the loop but occurs along the 
whole loop length, including the chromosphere. 

The considered model yields number of observational implications. The main 
feature is expansion of the column occupied by high-energy particles in radial di¬ 
rection. Taking into account recent progress in determination of geometric prop¬ 
erties of hard X-ray sources ( Kontar et ai, 20l0{ [Battaglia and Kontar, 2011 1, 
this feature can be used as an observational test for this sort of models. However, 
more simulations are needed to investigate this reconnection and acceleration 
scenario with wider range of parameters. 
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